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General properties of the tunnelling-charging Hamiltonian 
of a Cooper pair pump are examined with emphasis on the 
symmetries of the model. An efficient block-diagonalisation 
scheme and a compatible Fourier expansion of the eigenstates 
is constructed and applied in order to gather information on 
important observables. Systematics of the adiabatic pumping 
with respect to all of the model parameters are obtained and 
the link to the geometrical Berry's phase is identified. 

Adiabatic transport of single electrons in arrays of small 
metallic tunael junctions has been widely studied in re- 
cent years.Eli3 In the Coulomb blockade regime phase- 
shifted gate voltages have been used to induce a dc cur- 
rent / = —nef, where n is the number of carried elec- 
trons and / is the gating frequency. Normal-state pumps 
transporting single electrons have reached accucacy that 
can be considered for metrological applications.H Pump- 
ing of Cooper pairs has gained interest|-£iue to new ideas 
in quantum measuring and computingJj'Q 

A quantitative theory of pumping Cooper pairs in 
gated arrays of Josephson junctions when the envirc 
ment has negligible impedance has been presentedl 
The leading order pumped current is / « — 2e/[l — 
flAf(£j)^~^ cos 0], where qn is a constant and ej := 
E.]/Ec is the coupling strength. Here i?j and Eq := 
(2e)^/2C are the Josephson coupling energy and the 
charging energy, respectively. According to recent cal- 
culations the cos 0-dependent inaccuracy should be ex- 
perimentally observable at least in a certain frequency 
range.Q In this article the model for adiabatic pumping 
from Refs. |^ and |^ is examined in a thorough manner, 
concentrating on the symmetries and general features of 
the simplified system, thus extending the treatment to 
strong couplings and long arrays. The overall behaviour 
of the pumped current is explained and systematised, but 
not rigorously proven. 

A schematic view of a Cooper pair pump and the ideal 
operation of the gate voltages Vg,k are shown in Fig. ||. 
On any of the TV legs of a cycle at most two gate voltages 
are changed. The tunnelling-charging Hamiltonian 



H = He 



H, 



(1) 



neglects the quasiparticle tunnelling and other degrees of 
freedom. The full set of model parameters are ej, total 
phase difference over the array </>, the relative junction ca- 
pacitances c, where Ck ■= Ck/C and J2k=i ^k^ ~ N/C, 
and the (normalised) gate charges q :— {qi, . . . ,qN-i}, 
where qk := — Cg^fcV^,fc/2e. For homogeneous arrays 
Ck ■— 1 and the inhomogeneity is quantified by Xjnh = 



lY^ki^k^ - ly/N]^/"^. The model Hamiltonian is diago- 
nal with respect to 0, which is assumed to be fixed, but 
could be controlled by an external bias voltage accord- 
ing to d(p/dL= —2eV/h and is subject to any voltage 
fluctuations.LiB The conjugate variable M, the average 
number of tunnelled Cooper pairs, is undetermined in 
the present model. 
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FIG. 1. (a) A schematic drawing of a gated Josephson array 
of A'^(> 3) junctions. The Ck and Cg,fc are the capacitances 
of the junctions and gates, respectively, b) The optimal oper- 
ation of the gate charges q. Each cycle carries approximately 
one Cooper pair through the array, when £j <g 1. 

The matrix elements of the charging Hamiltonian He 
are given by the capacitive charging energy 



{n\Hciq)\n)^ = Ec 



N 2 1 / ^ 

V^- — IV — 



(2) 



where the number of Cooper pairs on each island is given 
by n. The quantities Vk, k — 1, . . . , N , are a solution of 



Vk 



Vk+i = n.k - qk- 



(3) 



Tunnelling of one Cooper pair through the fcth junction 
changes \n) by 6k, where the non-zero components are (if 
applicable) {Sk)k = 1 and {dk)k-i = ^1- The tunnelling 
Hamiltonian then reads 



H, 



N 

E 

n,k—l 



CkEj 



(|n + 4)(n|e'^/^+ H.c). (4) 



The supercurrent fiowing through the array is determined 
by the supercurrent operator 



1 



Is 



-2e/n){dH/d(l>), 



(5) 



a Gateaux derivativeS of the full Hamiltonian. By chang- 
ing the gate voltages adiabatically along a closed path 
r, a charge transfer Qtot := Qs + Qp is induced. The 
pumped charge, Qp, depends only on the chosen path, 
while the charge due to direct supercurrent, Qs, also de- 
pends on how the gate voltages are operated. If the sys- 
tem remains in a adiabatically evolving state |m)n^Jie 
total transferred charge, Qtot, in units of — 2e reads 



dErait) 



dt~2 



E 



Im 



{m\h\l){l\dm) 



El — E„ 



(6) 



where \dm) is the change in |m) due to a differential 
change of the gate voltages dq. 

Intermediate states \l) and energy denominators can 
be removed by rewriting the off-diagonal matrix element 
between stationary states as 



{m\dH/dcj,\l) = [El - Ernm{m\)/dm. 



(7) 



and using the identity 10(^1 — 1- The canonical rep- 
resentation M = —id/d4) connects the pumped charge to 
the average number of tunnelled Cooper pairs by 



Qp = 2 ^ Re {m\M\dm 



(8) 



i|dm).y The connection dQ-p{(l)) - 



A single eigenstate for phase differences (j) and 4> + d4) 
is required for each integration point. This expression 
identifies M as the link between Qp apd the geometrical 
Berry's phase, 7m(r) = 

was given without an explicit identification in Ref. ^ and 
was mentioned in Ref. For ej — > only two charge are 
of importance during each leg. From the wave function 
|m) = [(l-a2)i/2^ ^gj^/wjT^ vfhere a : ^ 1, one obtains 

{m\M\dm) = d{a'^)/2N and Qp = 1 for a full cycle. 

A reference state |no) induces a convenient labelling 
of the charge states. Each state is denoted by in- 
tegers {yk}k=i^ < Yfi := Y^kVk < such that 
n = fio + J^kVk^k- The numbers yk tabulate the num- 
ber and direction of tunnellings from uq to n. The dis- 
tance between charge states is defined by d{ni,n2) 

min (Er=i IVk^ - 2/f + ?| : ? e Z) . If 7x2) = 1 

(= I), then |ni) and \n2) are (Ith) nearest neighbours. 

A change of basis t7{|n)} = {|e*'^^"/^n)} yields a 
Hamiltonian matrix H = UHU~^ with properties H{(f) + 
27r) = H{(j)) and H{-(j)) = {&{(/)))*. Thus eigenstates 
and eigenvalues of H are periodic under 2tt, or state la- 
bels change cyclically. Usually the former happens, so 
the supercurrent in a stationary state is described by a 
Fourier sine series 



S/(N,E.,,<t,,c,q) '■- 



YaLi^i su\{l(j)). 



(9) 



The ground state supercurrent behaves differently only 
at the so-called resonance points, where the ground state 



becomes degenerate for = tt -I- 21-k, I G Z. For ho- 
mogeneous arrays resonance points are located at g = 
n±(l,...,l)/7V, where n is arbitrary. The correspond- 
ing poles in the Berry's phase give raise to Qp, and ej 
determines which poles are important. 

Due to symmetry of the representation H the original 
amplitudes in 27r-periodic states are given by 



= E 4, 



(10) 



where real Fourier coefficients at ^ are fixed by the gauge 

condition a'tf :— |a^',"^|e*'^-^"'/^ for the charge state \n'). 
The averaged number of tunnelled Cooper pairs. 



M:=j:nAl + yr-t/NKal,) 



(11) 



is unique up to a gauge-dependent integer. The apparent 
contradiction between sharp phase difference combined 
with sharp value of M is an artefact due to gauge-fixing. 
Stronger coupling increases the variance of M . Disconti- 
nuities in gauges certainly occur on any closed path en- 
circling an odd number of resonance points. A gauge is 
unstable near a discontinuity, but everywhere away from 
resonance points many valid gauges exist. Especially on 
the gating path depicted in Fig. 1, the dominant charge 
states on each leg give very stable gauges. 

The pumped charge, Qp, can be evaluated using a 
gauge-independent differential expression 



E E 

l'=0 n,l=-c 



2{l + YrjN) 
1 + 5vo 



I' {an,idafi. 



i+v 



an,i+i'dan,i)]cos{l'(f)), (12) 



where the (/)-independent average is simply dM. Due to 
the normalisation of Im)^, the coefficients are orthonor- 
mal, i.e. J^n i o,n,iO'n,i+i' = Si'o, which cancels the terms 
multiplying the full differential by T. Expression (^2|) in- 
dicates that the averaged charge transfer for a full cycle 
is is exactly — 2e, regardless of the inhomogeneity of the 
array or reasonable deformations of the gating path. 

The tunnelling-charging Hamiltonian can often be 
block diagonalised with the following transformation. 
Let an orthonormal basis {\s)} span the Hilbert space TC 
and the matrix elements of a Hamiltonian H be hgs' ■= 
{s\H\s'). Choose projection operators {Pi} by Pi — 
^di(<oo) require that J^ij PiPj = Et Pi 1- 

If all of the row sums 



W,, 



ij,k 



■ 1 ^ikik' 



fc = 1, 



(13) 



are independent of k, the Hamiltonian commutes with 
the projection operator P :— J2i where {ipi) = 

EfeLi l*fe)- Thus H can be written as a direct sum 
H = Hp ® Hi^p with matrix elements of Hp given as 

mmi) = W.,{d,/d,Y/^ = WUd./d^f^. (14) 
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The block-diagonalising transformation is not specific to 
the present problem, as it amounts to extracting states 
with a specific symmetry from all of the system's eigen- 
states. The ground state of He + Hj is always an eigen- 
state of Hp, because ai'^~'^ > for all q and n. The 
Fourier expansion ( p!o| ) can be used simultaneously if each 
subspace label, i, corresponds to a single value of Y,i. 

For homogeneous arrays at q = Hq all junctions 
are indistinguishable in terms of the charging energy. 
The subspace labels are of the form rhg with dm^ — 
Nl/{Ylj'^^ rrijl). This subspace contains charge eigen- 
states carrying the label 

,-(z(^...,z^),4^\...,zr-),...,z(2r-)), (15) 

where Zi < ■ ■ ■ < Zj_^^^ and mi + • • • + mj^,^^ = N, or any 
distinct permutation of y. Thus only the number and 
multiplicity of tunnellings is of importance. 

On the gating path the gate charges can be written 
as (f (1 — x)no + xfii, where |no) {\ni)) is the initial 
(final) optimal charge state for the given leg. Thus one 
junction becomes distinguishable and subspace labels can 
be chosen as zq; to^, where fhz refers to the remaining A'^— 
1 components, < zg + X^ilii^ < ^ and (izo;mj = 

(A^-i)!/(n^=rm.!)- 

The diagonal and non-zero off-diagonal matrix ele- 
ments of Hp are the common charging energies and 

- (ejm,/2)e±**/^(d(,„;)^,/d(,,0,y.;,)'/' (16) 

with mo = 1 when applicable, respectively. The block- 
diagonalised matrices are sparser than the original ma- 
trices and the amplitudes in eigenvectors are multiplied 
by (c?(zo;)m5)^^^ due to combining of several amplitudes 
into one. A further block diagonalisation is possible at 
q = nQ and q — (no + ni)/2 if is a multiple of tt. 

A truncation picks the charge states required for reli- 
able evaluation of eigenstates and observables. First an 
initial truncation, a set of states Bi = {\nj)}, is chosen. 

Its extensions, bases contain all neighbours up to 

and including Ith nearest neighbours of each and every 
\nj). For non- ideal cycles and/or inhomogcncous arrays 
the initial truncation is the "b basis" of Ref. ^ which re- 
produces leading order supercurrent and inaccuracy. A 
"c basis" is just the first order extension of a "b basis" , 
i.e. "bbasis(i)". 

For ideal cycles and homogeneous arrays the optimal 
Bq truncation corresponds to 3A^ — 2 labels, that is 
0; {j, Af-l-j)(o.i), "1; (j, iV~l-j)(_i,o)", J = 0, . . . , TV- 
1, and 1; (j,iV - 1 - j)(o,i), j = 1, . . . , TV - 2. These 
3.2^-1—2 charge states actually contribute to the leading 
order supercurrent and inaccuracy. More restricted, few- 
state truncations Bi = {|no)} and B2 = {|?^o)j l^i)} are 
of use especially when the ground state energy is sought. 

The efficiency of the block diagonalisation for bases Sq''* 
is shown in Table ^ where the number of labels is com- 
pared against the number of included charge eigenstates. 
Examples of bases B^ ^ and iJj'' are also given. 



TABLE I. The number of subspace labels (si) and included 
charge states (ch) for bases Bq'' with several A'' and examples 
of bases B['^ and 



^0/1/2 


si. 


ch 


si, 


ch 


4o 


522 


182 


1872i4 


9572 


5o 


822 


626 


3082ii 


46400 


60 


1142 


1918 


34689 


1.56- 10^ 


7o 


1462 


5428 


39418 


5.22 • 10^ 


80 


1782 


14498 


36487 


1.30 • 10*^ 


9o 


2IO2 


37082 


28996 


2.44 • 10*^ 


81 


695 


12331 


1472719 


6.15 • 10^ 


72 


2485 


9452 


12842i4 


2.19 • 10^ 



As a concrete example, take basis i?} and the sub- 
space labels iVo, (iV — 1, l)(o,i) and (1, TV — 1) (0,1), stand- 
ing for 2TV -I- 1 charge eigenstates at q — ftf). The block- 
diagonalised Hamiltonian obtained from Eqs. (^4|) and 
(|l6| ) is given by 

/ K*^ Kn \ 

H^A = Km {N - 1)/TV 2Sm3K*^/VN (17) 
\K*M 26n3Kn/VN (TV-1)/TV / 

where Kn := -(ej/2)e"^/^\/TV. For TV = 3 the lead- 
ing order supercurrent for ej <C 1 is already reproduced. 
Note that labels (TV — 1, l)(o,i) and (1, TV — l)(o,i) can be 
joined, if is a multiple of tt. 

The discussion has now lead to the main results of 
the paper, the general systematics of the pumped charge. 
For ideal gating sequences dAi = 1/TV for each leg due 
to symmetry. Corresponding integrated pumped charge 
reads 

Qp(TV, ej, c, leg) := 1/TV -I- ^Mcg cos(?0), (18) 

which is always positive. Replace 1/TV by dA4 for non- 
ideal cycles. When performing numerical calculations, 
the basis must contain all of the important charge states 
and the number of angles used in fast Fourier transform, 
2^, must be large enough not to allow misidentification 
of non-negligible components aaj and a^j j_|_2j as a single 
component. Failure to satisfy these requirements causes 
systematical error due to incorrect amplitudes and spuri- 
ous interference between Fourier components in Eq. (^) , 
respectively. The integration points should be chosen so 
that the magnitude of norms || 1^™)]! is neither too large 
nor varies too much. 

The rate of convergence for each cos(/'(/))-dependent 
term in Eq. ( p^ behaves as l/(#steps)^. For a reason- 
able choice of integration points 300-500 steps per leg, 
much less than used in Refs. 0and|^, usually gives relative 
precision of the order of 10 (disregarding the system- 
atical error) for coefficients with magnitude greater than 
10^^. This greatly diminishes convergence problems due 
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to smallness of dari^i for large bases and large number of 
angles. The expression 21' afi^idan^i+i' has been used as 
the latter part of Eq. ( |l2| ) for many of the data points. 
In this case the rate of convergence behaves as l/#steps. 
Bases -B2''' have also been used when calculating Qp. 

TABLE II. The first four ratios Pi and the limiting value 
P for £j = 0.4 as function of A'' and for A'' = 4 as function of 
ej. Subscript of /3 in the last column indicates the coefficient 
Pi used when estimating the iimiting value. 



(ej)iv 


Pi 




/?3 


Pi 


{P)i 


0.43 


0.8802 


0.7782 


0.7470 


0.7347 


0.72260 


0.44 


0.7250 


0.5706 


0.5354 


0.5225 


0.50925 


0.45 


0.5494 


0.3842 


0.3551 


0.3451 


0.33520 


0.46 


0.3844 


0.2393 


0.2188 


0.2121 


O.2O615 


0.47 


0.2516 


0.1401 


0.1272 


0.1231 


O.I2O11 


0.48 


0.1563 


0.0785 


0.0708 


0.0685 


0.0679 


0.49 


0.0934 


0.0427 


0.0383 


0.0370 


0.0367 


0.4io 


0.0542 


0.0227 


0.0203 


0.0196 


0.0197 


0.044 


0.0188 


0.0096 


0.0095 


0.0095 


0.00954 


0.24 


0.3444 


0.2164 


0.2049 


0.2025 


O.2OI14 


0.64 


0.8853 


0.7833 


0.7469 


0.7296 


O.7OI40 


0.84 


0.9475 


0.8876 


0.8598 


0.8444 


0.80880 


I.O4 


0.9738 


0.9390 


0.9198 


0.9079 


0.870i2o 


1.24 


0.9860 


0.9654 


0.9526 


0.9441 


O.9O8160 



ferred charge Qp = X^icg Qp.icg for inhomogeneous ar- 
rays does not have to satisfy /?; > /3/+i, although always 
—bi+i/bi < 1. This symmetry is also broken by non-ideal 
gating sequences, even for single legs. 

In order to conclude, the properties of the tunnelling- 
charging Hamiltonian of a Cooper pair pump have been 
examined using an efficient block-diagonalisation scheme 
and a compatible Fourier expansion of the eigenstates. 
Explicit enforcement of the model symmetries produces 
strong systematics of the pumped charge, even if the 
structure of the Fourier coefficients {bi,\cg}'i^o ^'-'^ 
haustively proven. These properties are possibly related 
to the dynamical algebra of single and coupled Josephson 
junctions described in Ref. |l^, which offers a complemen- 
tary view of the present problem in case of a supercon- 
ducting loop. 
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The sequence of the Fourier coefficients {6i,icg}^0' 
where bo — 'i./N , is alternating and decreasing in magni- 
tude. Furthermore, the ratios Pi :— |5//6;_i|, Z = 1, 2, . . . 
form a decreasing—sequence with limiting value [3 := 
limj^oo /9; > 6o/3i-li3 Because dQp/dcj) > in the range 
4> G (0,7r), Qp is bounded from above by [2(1 — (3i)~^ — 
1]/N and from below by 1/A^ - 2f3i/{N + 2Pi). For fi- 
nite values of £j the ratios /?/ and (3/ (3i are monotonously 
increasing functions of ej with limiting value of unity, be- 
cause Qp{(t> = 0, £j — > 00) \ 0, which enforces a stricter 
limit for the ratio (3/ [3i. This behaviour is explained by 
increasing long-range (high-Z') correlations in Eq. (|l^) or 
actually in the state itself. On the other hand, the ratios 
are monotonously decreasing functions of N for fixed ej 
as correlations are weakened in longer arrays. All of the 
above-mentioned features are clearly shown in Table ||, 
where (3i for the homogeneous case are depicted as funjG-. 
tion of A'^ at ej = 0.4 and as function of ej for N = 4.t2l 
These properties of the pumped charge are rather robust 
against relatively small systematical errors . 

For small values of £j <C 1 one has /3;>i ~ 6o/?i and 
I3i ^ {Nej/2)^^'^[N{N -1)/2{N -2^]. The effects due 
to inhomogeneity, that is the relative sizes of 6i,iog, can 
estimated by fixing the the leg index r in the inhomo- 
geneity prediction, Eq. (37) of Ref. ^. The correspond- 
ing results work quite well, only slightly overestimating 
the ratio between largest and and smallest coefficients, 
even for large inhomogeneities Xjnh- The total trans- 
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